arXiv:1507.06177vl [cond-mat.dis-nn] 22Jul2015 


cpi uiaiL 


Correspondence between many-particle excitations and the entan¬ 
glement spectrum of disordered ballistic one-dimensional systems 

Shaul Leiman^, Ariel Eisenbach^ and Richard Berkovits^ 

^ Department of Physics, Jack and Pearl Resnick Institute, Bar-Ilan University, Ramat-Gan 52900, Israel 


PACS 61.43.-j - Disordered Solids 
PACS 73.21. Hb - Quantum Wires 

PACS 03.65.Ud - Entanglement and Quantum Non-locality 

Abstract - Using exact diagonalization for non-interacting systems and density matrix renormal¬ 
ization group for interacting systems we show that Li and Haldane’s conjecture on the correspon¬ 
dence between the low-lying many-particle excitation spectrum and the entanglement spectrum 
holds for disordered ballistic one-dimensional many-particle systems. In order to demonstrate the 
correspondence we develope a computational efficient way to calculate the ES of low-excitation of 
non-interacting systems. We observe and explain the presence of an unexpected shell structure 
in the excitation structure. The low-lying shell are robust and survive even for strong electron- 
electron interactions. 


Introduction. — Concepts developed in the field of 
quantum information are gaining a foothold in condensed 
mater physics m- One of the most influential concepts 
is related to quantifying the entanglement between two 
parts of a system. Specifically, a many particle system in 
a pure state, may be divided into two regions A and B. 
The entanglement between the regions A and B can be 
quantified by different measures such as the entanglement 
entropy, Renyi entropy and entanglement spectrum (ES), 
connected to the reduced density matrix of area A, pA or 
B, pb- 

The canonical method for utilizing the information em¬ 
bedded in the eigenvalues Xf of pA is the entanglement 
entropy: Sa = ~ ■ The ES is constructed 

out of the set of eigenvalues {Af} by the transformation: 
{ef = — InAf}. Several years ago Li and Haldane [3] 
came up with an intriguing conjecture regarding the ES 
of a many-body state. They show that the ES of a par¬ 
titioned fractional quantum Hall p = 5/2 state resembled 
the minimal edge excitation spectrum, and suggest a con¬ 
nection between the properties of the ES and the topolog¬ 
ical order of this state. 

This idea led to many additional studies EMI], mainly 
focused on systems with topological behavior, suggesting 
that the low-energy ES distribution shows some correspon¬ 
dence to the true many-body excitations (MBE) of the 
partitioned segment (region A). Since the reduced density 
matrix of a region is adiabatically connected to the MBE 


of the disconnected region, one may expect that the exact 
ground state eigenfunction of the whole system encodes 
information on the sub-system’s low lying excitations. 

In this letter we would use Li and Haldanes’ conjecture 
in order to investigate properties of the excitation spec¬ 
trum of fermionic disordered systems. The effect of disor¬ 
der on the ES of edge states in topological insulators has 
been previously studied . Properties of several low 

lying ES levels have been used to identify a metallic phase 
in one-dimensional Bose-Hubbard models [23]. Neverthe¬ 
less, a direct comparison between the properties of many 
particle excitations in a weakly disordered systems and the 
ES of a corresponding segment is still lacking. 

The field of MBE statistics has a long history, and con¬ 
tinues to draw interest in diverse areas such as cold atoms 
in the presence of quasi-periodic potentials [25H28] . and 
the many body localization (MBL) transition [29l - l3T| . For 
disordered systems it is natural to seek knowledge on the 
statistical properties of spectrum. Specifically, we are in¬ 
terested in answering the question: Does the ES exhibit 
statistical properties corresponding to the MBE spectrum. 
For the disordered single particle spectrum there is an ex¬ 
tensive literature on the statistical properties of the energy 
spectrum in the localized, critical, diffusive and chaotic 
regimes. Different energy spectrum and wave function 
statistics, depending on whether the disordered system has 
time reversal symmetry (Gaussian orthogonal ensemble 
(GOE)), broken time reversal symmetry (Gaussian uni- 
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tary statistics (GUE)), spin-orbit interactions (Gaussian 
symplectic statistics (GSE)) The energy spec¬ 

trum statistics can be used to identify the Anderson local¬ 
ization transition [35) . 

The statistics of MBE in disordered interacting sys¬ 
tems have an interesting twist. Eor non-interacting many- 
particle systems the level spacing MBE distribution is ex¬ 
pected to follow the Poisson distribution for excitation en¬ 
ergies above the second spacing, without depending on the 
single-level distribution [3B]. On the other hand, once re¬ 
pulsive interactions between the particles are considered, a 
transition to the Wigner distribution for higher excitations 
is observed [37H4T] . The main difficulty in studying this 
transition is that exact diagonalization needed to study 
excited states is limited to very small systems. Here, the 
Li and Haldanes’ conjecture can come to the rescue, since 
as we shall demonstrate below, one may extract the ES 
up to a few hundreds of states. Thus, if indeed there is a 
connection between the ES and The MBE, the low level 
excitations of a rather large many-body systems are nu¬ 
merically available. 

In this letter we shall explicitly demonstrate the corre¬ 
spondence between the excitation spectrum and the ES, 
for non-interacting one-dimensional disordered systems in 
the chaotic regime. In order to facilitate the calculation of 
the ES for large systems, we present a numerical method 
based on the correlation matrix eigenvalues, which could 
be used for non interacting systems in any dimension. We 
show that the average entanglement spectrum level spac- 
ings (ESLS) show a shell structure, with a large average 
spacing appearing according to the combinatorical parti¬ 
tion function pirn). The ESLS follows the Poisson statis¬ 
tics for the small spacings, while it shows a narrow distri¬ 
bution for the large spacings. Adding interactions removes 
the shell structure for strong interactions at higher levels, 
but the shell structure is rather robust for the low-lying 
portion of the ES, corresponding to the MBE which are 
close to the Eermi energy. 

Model. — As an example of a disordered many- 
particle system, we consider a spinless ID electrons system 
of size L with repulsive nearest-neighbor interactions and 
on-site disordered potential. The Hamiltonian is given by: 

L L-l 

^ ^ -b h.c.) (1) 

i=i 

+ UY^[c]cj--) (c]+iCj+i - -), 

i=i 

where is the on-site energy drawn from a uniform dis¬ 
tribution [—IE/2, IT/2], c] is the creation operator of an 
electron at site j, and f = 1 is the hopping matrix element. 
The repulsive interaction strength is depicted by 1/ > 0, 
and a positive background is considered. 

Non interacting electrons. — All single-electron 
states of such a ID system are localized [13], with states 



Fig. 1: The single electron level spacing distribution P{si) for 
states i = 173 ... 179 around the middle of the band for a 
L = 350, W — 0.3 wire. For this strength of disorder the 
system is disordered chaotic L < G s-nd the distribution is 
peaked around the mean (where the level at the middle of the 
band i = 176 is more peaked than the others). The Poisson 
and GOE distributions are provided for comparison. 


at the middle of the band having a localization length 
^ « 105/IE^ jUj. We consider a case where the disorder 
was chosen as IE = 0.3 and L = 350, i.e., ^ ^ 1100 > L. 
For this case the single electron energies an eigen¬ 
states tpi are readily available via exact diagonalization. 
Unless otherwise specified, we perform our analysis over 
10, 000 realizations of disorder. First we consider the sin¬ 
gle electron level-spacing distribution, for the z-th spac¬ 
ing 5i = Ei+i — £i. The unfolded spacing is defined as 
Si = Si/{Si), where (...) denotes averaging over the differ¬ 
ent disorder realization. The single electron level-spacing 
will follow the Poisson statistics (Ppoisson(si) = exp(—s^)) 
as long as L 3> G while in the metallic (diffusive) regime it 
should follow the Wigner (GOE) distribution (PGOE('Si) = 
{TTSi/2)exp{—TTSi /4:)) [35]. For a one-dimensional disor¬ 
dered system there is no genuine metallic regime since 
L < ^ ^ i (where i is the mean free path) and there¬ 
fore the system crosses from a localized to a disordered 
ballistic regime, resulting in a narrower distribution con¬ 
centrated around the mean {si = 1). This can be clearly 
seen in Fig. [TJ where the distribution of the single electron 
level spacings close to the middle of the band of a ballistic 
system is presented. 

The MBE of the system can not be obtained by exact di- 
agonaliztion since the size of the many-body Hilbert space 

grows as ^ N ) ’ half-filling {N = Ll2 = 175), 

where N is the number of particles, is of order 10^°^. Thus, 
another tack is needed. For the non-interacting case one 
may calculate the low-lying excitations by considering the 
different permutations of the single electron occupations, 
TT-i = 0,1, of the single electron states. Here, each permu¬ 
tation which obeys N = rii is a many-body state, with 
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Fig. 2: The properties of the many-body excitation spectrum 
for 10000 realizations of the non-interacting system {U = 0) 
whose single electron level spacings is depicted in Fig. [T] Top 
panel: Level spacing distributions for the eight lowest many- 
particle levels. Lower panels: Average spacings as function of 
the level number. The vertical line depict the positions of the 
peaks predicted by the partition function p{n). 


the total many body energy E = These many- 

body states may be enumerates by arranging the energies, 
such that El < E 2 < E 3 < .... Of course there are still an 
astronomical number of these states and covering them all 
is impossible. Nevertheless, if one constrains the search to 
low lying excitations, one can limit the the permutations 
to p (of order p ~ E/5) single electron states around the 

Fermi energy. Thus, only ^ ^ MBE are considered. 

Similar to the single particle case the MBE level spacing 
Ai = Ei+i — Ei can be extracted, as well as the average 
(Ai) and distribution P{Si), where Si = Ai/{Ai). 

The results for the average MBE spacings and their 
distribution are plotted on the right panels of Fig. [51 
It is immediate apparent that the average MBE spac¬ 
ing exhibits a shell-like (also known as magic num¬ 
bers in context of nuclear physics) structure. Spacings 
1,2,4,7,12,19,30,45,67... are significantly larger than 
their neighboring spacings. The reason for this behav¬ 
ior is combinatorical. The large spacing corresponds to 
an excitation different from the previous one by a sin¬ 
gle electron moving to the next available single-electron 
state. On the other hand the typical smaller spacings 
corresponds to the movement of several electrons be¬ 
tween different single-electron states. For equal single- 
level spacings (i.e., 5i = <5), the MBE can acquire ener¬ 
gies of md above the ground state where m = 1 , 2 ,.... 
The degeneracy of the states with energy mS depends 
on the the number of ways of writing the integer m as 
a sum of positive integers, where the order of addends 
is not significant. This corresponds to a problem in 
combinatorics knows as the partition problem, and the 
answer is given by the partition function p(jn) (which 



Fig. 3: The properties of the many body entanglement spec¬ 
trum for an La ~ 350 region out of a L = 700 system, with 
U — 0, averaged over 10000 realizations of disorder. Top panel: 
Entanglement level spacing distributions for the eight lowest 
levels for Na ~ La/'^- Lower panel: Average spacings as func¬ 
tion of the entanglement level number. Blue circles correspond 
to Na = LaI‘2,. Red squares to an additional average over 
LaI‘^ — 6 < Na 'C La/^^ 6. 


is not the thermodynamical partition function), which 
has a recurrence relation pirn) = ~ 

ki3k — l)/2) + pirn — fc(3fc + l)/2)] [45]. Thus, p{m) for 
m = 1,2,... are 1,2,3,5,7,11,15,22,30,42,.... Hardy 
and Ramanujan have devised an asymptotic form pirn) ~ 
(4m-\/3)“^ exp(7rY^2m/3) [35]■ In the case where the sin¬ 
gle level spacing is not constant the degeneracy is lifted, 
but as long as the variation is not too large, the shell struc¬ 
ture remains and the MBE level spacing between a closed 
shell and the next open one is relatively large, as can be 
seen in Fig. |5| Thus, the number of MBE states between 
the j-th MBE shell (i.e., large level spacing) and j + l shell 
is p(j), in agreement with the numerical results. 

The MBE spacing distributions T’(S'i) also retain the 
shell structure, where the distributions corresponding to 
the large spacings are similar to the single electron spac¬ 
ings (Fig. |T|), i.e., narrow and peaked around the mean, 
while for the other spacings follow the Poisson distribu¬ 
tion, expected when the non-interacting MBE differ by 
more than a single electron occupation [36l|37] (see Fig. 

E]). 

Does the ES of a hnite segment of a disordered system 
exhibit similar behavior as might be expected from the Li 
and Haldanes’ conjecture [5]? First we have to calculate 
the ES, i.e., the eigenvalues of the RDM. The density ma¬ 
trix renormalization group (DMRG) |47l|48] is a natural 
candidate for calculating the ground state of disordered 
interacting ID system and the corresponding eigenvalues 
of the reduced density matrix. We will use it for the in¬ 
teracting case. For non-interacting systems, one would 
expect that it is possible to extract the eigenvalues with¬ 
out to resort to DMRG. Direct calculation of the RDM is 
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impossible since its size is of order 2^-^ (where La is the 
size of region A), so a different approach is needed. 

Reduced density matrix eigenvalues using the 
correlation matrix. — In order to calculate the ESLS 
for the non-interacting case we use the connection between 
the correlation matrix (CM) and the RDM [?5] . The CM 
between two sites in region A is a unitary matrix given by: 

C'mn = [Tr{pAc\nCn}]^ = T T {clcm P a}^ C nm ■ (2) 

One may diagonalize the CM and write the trace in terms 
of the eigenvectors. Thus, 

Tr{pAa\ap} = Vq5qp (3) 

where aj is the creation operator of the g-th eigenvector 
and Vq is the corresponding eigenvalue. The RDM for non¬ 
interacting electrons can be written as an product state of 
local density matrices [32]: 

PA = Pi® P2® PLa (4) 

where Pq is the density matrix associated with the g’th 
eigenvalue of the CM. Moreover, due to the fermionic na¬ 
ture of the particles, pq has to be diagonal. 

The relation between Vq and pq is easily extracted [42] : 


Aqq = Tr{pAa\aq} = Tr{pqa\aq} = 



Since pq are independent any eigenvalue of the RDM 
can be constructed by multiplying a permutation of ei¬ 
ther (1 — Vq) (no particle) or Vq (one particle) occupy¬ 
ing the g-th state. For the case where no particles oc¬ 
cupy region A {Na = 0), the single eigenvalue of the 
RDM is = ng=i(l - '^q)- For the Na = 1 

the RDM has La eigenvalues, where the i-th eigenvalue, 
= (1 - vi){l - V 2 ) ■ ■ .Vi.. .{I - vla)- Generalizing 
to any Na = P particle block in the RDM the first eigen¬ 
value is constructed hy = vi... Vp{l — Vp+i)(\ — v^^), 

while the other ^ ^ ^ permutations define the rest of 

the eigenvalues. Since usually many of the CM eigenvalues 
are either extremely small or very close to one, it is possi¬ 
ble to reach an accurate estimation of A, with significantly 
fewer permutations. 

The ES obtained for a wire of length L = 700 and 
La = 350 occupied by Na particles and averaged over 
10,000 realizations of disorder, are presented in Fig. |31 
The occupation is Na = LaI‘^ or averaged for values of 
LaI‘2. — 6 < Na < LaI^^. 6. This size was chosen in 
order to correspond to the sizes of the system for which 
the MBE were calculated. For the ES we perform the 
transformation = — ln(Af^'^) and calculate the ESLS 
, the average spacing (A^) and distri¬ 
bution P{Si), where Si = Ai/{Ai). Indeed, the low-lying 


ESLS average as well as the spacing distributions show 
the same general features shown by the MBE. The shell 
structure is reproduced for spacings 1, 2,4, 7,12, small de¬ 
viations in the peak positions appear for the higher spac¬ 
ings at 19 and 30, and is almost wiped up at peaks higher 
than 45. The general behavior of the ESLS follows the 
distribution exhibited by the excitation spacings with nar¬ 
row distributions for spacings 1, 2,4 and 7, while the other 
spacings follow the Poisson distribution. Nevertheless, de¬ 
tails such as the exact width of the distribution vary be¬ 
tween the excitation and entanglement spacings. 

The influence of interactions. Once interactions 
are added, our previous arguments do not hold anymore. 
Nevertheless, following the ideas leading to the Fermi liq¬ 
uid picture and quasi particles, one expects that for not 
too strong interactions and close to the Fermi energy the 
non-interacting shell picture will continue to provide a 
good description of the system. 

As we can not use the combinatoric approach to calcu¬ 
late the interacting MBE energy levels, and neither exact 
diagonaliztion nor DMRG can provide more than a few ex¬ 
cited states, we will use the correspondence between the 
MBE and the ES of a finite section shown for the non¬ 
interacting case to study the evolution of the spectrum as 
function of electron-electron interaction strength. Thus we 
turn on the interaction U in Eq. m and use DMRG to 
calculate the ES over 100 realizations of the same disorder, 
wire length, size of region A, and number of particles as 
for the non-interacting case. The average ESLS A^, and 
distribution P{Si), are calculated and presented in Fig. 
131 For weak interactions U = 0.3, the shell structure of 
the average spacing is hardly affected. As interactions in¬ 
crease, the higher shells are washed out, until for U = 2.4 
only spacings 1, 2 and 4 remain larger than their neighbors. 
The entanglement spacing distribution is also transformed 
as function of interaction strength. For weak interactions 
{U = 0.3) the distribution P{S 3 ) is close to Poisson, while 
the distribution for the 4-th (large) spacings, P{S 4 ), is 
peaked around the mean, as for the non-interacting case 
(see inset Fig. SI). This distinction is blurred as the inter¬ 
actions become stronger and P{Si) for any i approaches 
the GOE distribution. This is in agreement with the ob¬ 
servations for the distribution of interacting MBE [5SMT| 
which show a transition to GOE statistics as interactions 
become stronger. Thus, it seems that the Li and Haldane 
conjecture holds even for the ES of interacting systems. 

Discussion. — Thus, Li and Haldanes’ conjecture on 
the correspondence between the low-lying MBE spectrum, 
and its ES is confirmed for non-interacting disordered sys¬ 
tems. Turning on interactions gradually change the ES of 
the system in a manner that fits the expectations for the 
behavior of the MBE. We conclude that ES is a promis¬ 
ing route to study MBE properties beyond the first few 
excitations, which are impossible to study by any other 
way for large systems. This could be a promising way 
to study phenomena occurring only for the MBE, such 
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Fig. 4: Averaged entanglement spectrum level spacings A„ as 
function of the level number n for different interaction strength 
U. Inset: The spacing, S, distribution for the 3rd (continuous 
line, empty symbols) and 4th (dashed line, full symbols) for 
U = 0.3,1.2 and 2.4. The Poisson (continuous black line) and 
GOE (dashed black line) are drawn for comparison. 

as the many-body localization transition [29]. Although 
still limited to low-lying excitations, the fact that higher 
many-body excitation may be probed improves the possi¬ 
bility to glean useful information on the transition. Indeed 
in the behavior of the shell structure as function of inter¬ 
action strength presented in Fig. [d] where as interaction 
increases lower shell peaks disappear, may show a glimpse 
of this phenomena. 

For non-interacting systems we have presented a com¬ 
putational effective method to extract the low-lying ES 
from the CM. This method could be used for the study 
of the ES for higher dimension, which is otherwise quite 
daunting. 

We have also shown a new shell (magic number) struc¬ 
ture appearing for disordered ballistic one dimensional sys¬ 
tems and explained its origin. The low-lying shells are 
robust and are not wiped out by interactions. It is inter¬ 
esting to understand whether this behavior survives for 
stronger disorder where the localization length is smaller 
than the length of region A. On one hand, since the disor¬ 
der is larger it is expected to wipe out the shell structure, 
on the other hand, the area in region A entangled with the 
rest of the system remains of order of ^ and thus shrinks 
with the disorder, i.e., one effectively samples a smaller 
sample. Also the study of higher dimensions may prove 
interesting. 

* * * 
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